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Abstract. We introduce a model which incorporate the subdiffusive dynamics 
and the ratchet effect. Using a subordination ideology, we show that the resulting 
directed transport is sublinear, {x(t)} ~ Jt^ , /3 < 1. The proposed model may be 
relevant to a phenomenon of saltatory microbiological motility. 

Thermal fluctuations alone cannot create a steady transport in an unbiased system. Luck- 
ily, microbiological realm operates far from equilibrium [1] , where directed motion can appear 
under nonequilibrium conditions [5]. The corresponding ratchet effect has been proposed as 
a physical mechanism of a microbiological motility [3]. The nonequilibrium conditions might 
induce another intriguing peculiarity of microbiological transport; namely, the anomalous dif- 
fusion [4] . The quasi-random wandering at the molecular scale can be characterized by a mean 
square displacement (msd), a{t) = {x'^(t)) — (a;(t))^, which, in many cases, follows a power law, 
(T(t) ^ i", a ^ 1, rather than the linear time dependence for a Brownian particle |5l6j . a > 1 
corresponds to enchance or superdiffusion, while a < 1 to subdiffusion [4I7| . 

Directed current and diffusion are generated by the same trajectory, and characterized 
by the first and the second moments of the same process being therefore strongly conjugated. 
In the case of the normal diffusion and the superdiffusion, due to finiteness of all statistical 
moments, a directed transport, if any, corresponds to a linear grow of the mean displacement, 
{x{t)) ~ t, with the corresponding current J = hmt_^oo x{t)/t [5]. 

For the subdiffusion a situation is less obvious. In this regime the motion is made up of 
periods of sticking events separated by fast jumps to a new position (see, f.e., Ref. [6]). Formally, 
the broad power-law distribution of sticking times, 

7/;(t)^<-i-^ </?<!, (1) 

leads to the subdiffusion with an exponent a = /3 [7]. The saltatory molecular motor's trans- 
port [9], rapid bursts of directed movements interrupted by pauses of variable duration, has 
been tracked within a cell by using microinjected fluorescent beads [TU]- The abovementioned 
observations call for a study of a ratchet transport mediated by anomalous trapping events. 

In this paper we propose a simple model which naturally incorporates both mechanisms, 
subdiffusion and ratchet effect, thus bridging two research lines that were so far basically dis- 
connected from one another. Using subordinated ideology [TTj we show that the subdiffusion in 
a ratchet potential results in the sublinear directed transport, < x{t) >~ Jt^ . The subordinate 
formalism enables us to reformulate the problem within the circle map's theory and to derive 
necessary and sufficient conditions for the directed current appearance. 

We start with the model which describes a dynamics of the overdamped particle exposed 
to the shot-noise, 

where aj gives the length and the direction of the corresponding step taking place at the 
time instant tj. We assume that the length Oj depends locally on a periodic potential U{x), 
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U{x + L) — U{x), and a noise ^, such that 

a,=-U'ix{t,))+at,). (3) 

The model (t^JEl) can be treated as the overdamped hmit of the standard model [12 . The 
corresponding process, x{t), can be considered not as a function of time t, but rather as a 
function of the number of steps, n. The dynamics of the system ([2]) can be represented as the 
noised circle map \n\ . 



Xj + f{Xn) +^n, (4) 



where f{x) = —U'{x) stands for the acting force. The process x{t) is subordinated [TT] to the 
map (HI), such that the time is governed by the linear map, 

tn+l ^tn + Atn. (5) 

In addition, here we assume that (i) the dispersion is finite, {aj) < oo, and (ii) the time 
between steps, Atj — tj^i — tj, is a random stationary process with the probability density 
function (pdf) tplAt). If tl){l\t) has a finite first moment, (At) < cx), and aj is independent 
on X and has a symmetrical distribution (Gaussian, Poisonian, etc) then for the time scale 
t » (At) we get normal Gaussian diffusion, < x'^{t) >~ t, which can be described by an 
ordinary Langevin equation [11] . If ^(t) has a divergent first moment, which is the case of the 
distribution ([T]), then we deal with the subdiffusion, where the corresponding msd's exponent 
isa = /3 H. 

The role of the sticking time reduces to the fact that the actual number of steps made up 
to the time instant t fluctuates, so the operational time n is a random function of the physical 
time t. This function, however, is monotonously nondecaying with t and thus allows a causal 
ordering of the events. For the pdf p(x, t) one has [TT] 

p(a;,i) =^VF(x,n)x„(t), (6) 

n 

where W{x,n) is the pdf for the iterated process ([4]), and Xn{t) is the probability to make 
exactly n steps up to time t. 

The asymptotic transport is {x{t)) = J-{n{t)), where {n{t)) — X]^o "'Xn(^)i and the current 
value J follows from transport properties of the map (|3]), {x{n)) ~ Jn. By using the Laplace 
transform in the time domain, it can be shown that {n{s)) = ip/s^l — ?/;), where ?A(t) is the pdf 
for sticking time. For the Poissonian process, ijj^s) ~ i//{s + i^), one can easily get {n{t)) — vt. 
For the power-law pdf V'(0 follows that {n{s)) w t^^s^^^^ and, finally, {n{t)) « 7Ff(i+^)7 so 
that 

where r[x) is the Gamma function. 

The subordination approach allows us to separate transport properties, precisely the value 
of the generalized current J, which follows from the map from the sublinear asymptotic 
([7]), which is governed by the sticking time pdf V'(0- 

As an illustrative example we consider here the two-harmonics potential force, 

f{x) = El sin(27ra;) + E2 sin(47ra; + 6*), (8) 

which transforms the system (j4]) into a ratchet version of a climbing circle map T2]. Without 
loss of generality, we chose here the Gaussian noise ^ with the dispersion rf. 

Certain symmetries of the Eqs.([2l[8]) need to be broken in order to fulfill necessary conditions 
for a directed transport appearance |15| . Suppose that there is a transformation which leaves the 
equation ([2]) invariant and changes either the sign of cc, a; — » — x, or invert time, t — > — i (but not 
both operations simultaneously!). Such a transformation maps a given trajectory into another 
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Fig. 1. (a) The bifurcation diagram and (b) the current value J as functions of E2 for the map (I4I8|) 
(-Bi = 0.5). Black line corresponds to deterministic case, rj^ = 0, and the grey line corresponds to 
noised case, rj^ = 0.1. Note that the current value in the last case is scaled by the factor 40. 



one, with the opposite velocity. If at least one of such symmetries exists, then contributions of 
the trajectory and its symmetry-related counterpart will cancel each other and the asymptotic 
current J will be equal to zero [15j . 

In order to fulfill the necessary condition for the dc-current appearance, the following sym- 
metries of the potential force should be violated: 



,f{x)^-f{~x) (or U{x)^U{-x)), (9) 
f{x)^-f{x + L/2) (or U{x)^-U{x + L/2)), (10) 

which are the reflection- and the shift-symmetry following to Ref . |15 and termed as " symmetry" 
and " supersymmetry" in Ref. 16 . Both the symmetries are violated when E2 ^ and 9 ^ kn. 

The symmetry violation is the necessary condition for the current appearance. The current 
value is determined by microscopic dynamical mechanisms |17j . It is reasonable then to start 
the analysis of map's transport properties from the deterministic limit, ^ = 0. Fig.l shows 
the bifurcation diagram and the current value as functions of the second harmonic amplitude, 
i?2. There is an evident relationship between multiple current reversals and different kinds of 
bifurcations, which is a general property of underdamped deterministic ratchets [18]. Here, at 
the overdamped limit, this relationship can be explained qualitatively. The current value is 
equal to the average force, 

J={f{x))^ [\xf{x)P{x), (11) 
Jo 

where P{x) — X]J^-oo + ^L), n G Z, is the reduced invariant density for the map ([4]- 
[8|). Since any critical bifurcation, like a tangent bifurcation (transition from a limit cycle to a 
chaotic regime) |12j . always accompanied by drastic changes of the invariant density P{x), such 
a bifurcation leads to a "jump" of the current value. 

The addition of a noise changes the rectification dynamics, but still the system generates a 
non-zero current (gray curve in Fig. lb). However, the invariant density even at a weak noise 
limit converges to the Boltzmann pdf [19], P{x) = exp{—U {x) / a^) , so that the integral in 
rhs of Eq. (jlip goes to zero. Thus, the current decays rapidly with the increasing of the noise 
strength (Fig. lb). 
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Fig. 2. The evolution of the ensemble mean displacement {x{t)) versus t for the different values of 
subdiffusive exponent: a = 0.5 (•), a = 0.7 (■), and a — 1.2 (♦). Lines correspond to asymptotics 
from Eq. The parameters are N = 10^ , Ei = 0.5, E2 = 0.8, 9 = -2.247r, and rj^ = 0.02. 



We furthermore assume that a time between subsequent kicking events is a random station- 
ary process with the pdf mQ. The asymptotic displacement can be written as {x{t)) = J-{N{t)), 
where the current value J follows from transport properties of the map {x{n)) ~ Jn. 

Formally, we get 

{x{t)) ~ je. (12) 

We consider now a large ensemble of noninteracting ratchets. The dynamics of each particle 
in the operational time frame is governed by the same map (jHH]), but the physical time t is 
different for different particles. Nevertheless, the subordination formalism allows to make casual 
ordering of events [11]. Thus we can calculate the displacement {x{t)) by using an ensemble 
averaging. In Fig. 2 we shown the evolution of the mean displacement, {x{t)) , for different values 
of the waiting time exponent a. 

Spatial distributions for different times are shown on Fig. 3. For the case of normal dif- 
fusion (Fig. 3a), a > 1, the evolution follows an universal Gaussian scaling, which has been 
found for weakly underdamped ratchets [20 . This Galilei invariant Brownian process, p{x,t) ~ 
--^g( ^^^* ), where g{x) stands for the Gaussian pdf, is very different from the Galilei variant 

subdiffusive ratchet regime |7j. The corresponding pdf is asymmetric with respect to its cusp- 
like maximum which stays fixed at the origin, and the plume stretches more and more into 
transport direction (Fig. 3b). This behavior is reminiscent of a subdiffusive dispersive transport 
under a constant tilting force |21j . 

Summing up, we have considered the new model which yields anomalous ratchet dynam- 
ics. Transport properties, such as direction and value of the generalized current stem from the 
ratchet-like periodic potential. The anomalous character of a kinetics is governed by the waiting 
time pdf. The proposed model provides further contribution to the studies of the microbiological 
transport. It sets up the link between ratchets and a power-stroke approach to a microbiolog- 
ical transport, still existing dichotomy ^22]. The model can also be useful for an analysis of 
experimentally detected saltatory motility in a cell by using of microinjected fluorescent beads 

m- 



^ For the generation of the random variable At with the pdf ([1} we have used the random variable 
5 with the uniform distribution on the unit interval, [0, 1], and the transformation At — Atc£,~^^'^ . 
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Fig. 3. Spatial distributions for the ensemble of A'' = 10^ single ratchets for (a) a = 1.2 and (b) a = 0.5 
for different times t. The inset shows single particle's trajectory. Other parameters are as in Fig. 2. 
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